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1 Introduction and Summary 



To a large extend motivated and guided by the gauge theory/string theory correspon- 
dence of Maldacena [1-3], Baier et.al [4] and Bhattacharyya et.al [5] recently formulated 
second order relativistic viscous hydrodynamics of conformal fluids. It was found that 
Mueller-Israel- Stewart (MIS) theory [6,7], often used in hydrodynamic simulations of 
strongly coupled quark-gluon plasma (sQCD) produced in heavy ion collisions at RHIC, 
does not properly account for all second order viscous corrections. The framework of 
this new viscous hydrodynamics was applied to bulk physics analysis at RHIC in [8]. 

In this paper we continue the program of computing the transport coefficients of 
strongly coupled gauge theories from the dual string theory holographic description. 
The motivation to do so, even though we lack a tractable string theory dual to QCD, is 
the amazing regularity of transport coefficients for a large class of gauge theory plasma 
at strong 't Hooft coupling [9-16], and the fact that sQGP appears to exhibit near- 
conformal dynamics. Specifically, we concentrate on finite 't Hooft coupling corrections 
to two second order transport coefficients of a CFT plasma — a relaxation time rn and 
K. These transport coefficients are the siniplest to compute as they can be extracted 
from the equilibrium correlation function^j. While we perform the computations in 
the string theory dual to A/" = 4 SU{Nc) superconformal Yang-Mills (SYM) theory, 
following the arguments in [160, the results are expected to hold for any conformal 
gauge theory plasma in four dimensions, and should provide a reasonable estimate for 
sQGR 

Following [4, 5] , the stress energy tensor T^'^ of a viscous conformal fluid to second 
order in derivatives of the four-velocity u'^ takes form 

T^"" = eu^u" + P {g"" + u^'u") + W , (1.1) 

where e, P are the energy density and the pressure, and the dissipative part U^'^ is a 
sum of the flrst (Hi) and the second order (H2) terms: 

H^'^ = Hf {VaUfs; {v}) + nr {V^VaUf,; {v-m, Ai, A2, A3}) . (1.2) 

The four-velocity dependence of the dissipative part of the stress-energy tensor is 
uniquely flxed by conformal invariance up to a single phenomenological coefficient 



^The coefficient Ai is necessary to describe the Bjorken expansion [18] of plasma and its finite 

coupling corrections will be discussed elsewhere, 
more detailed discussion is given in [17]. 
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to the first order in the derivative expansion, and up to five additional coefficients to 
the second order in the derivative expansion. In the former case, a single coefficient 
is a shear viscosity rj, while the latter includes a relaxation time rn, three coefficients 
Aj (describing viscous terms bilinear in four- velocities), and a K-term (describing vis- 
cous hydrodynamics in curved backgrounds). In the linear regime, and in Minkowski 
space-time, i.e., setting Aj = k = 0, the new viscous hydrodynamics of [4,5] reduces to 
MIS theory [6,7]. A crucial observation of [4,5] was that the MIS regime of a strongly 
coupled four-dimensional conformal gauge theory plasma is simply inconsistent. In a 
specific exampl^ of A/" = 4 SYM, and at infinite 't Hooft coupling A = Qym^c 
it was found that [19,4,5] 



oo, 



rj 1 2 — In 2 r] 

s 47r ' ^ 2nT ' ttT 
rj rj 

Ai = — — , Ao — ~ , A^ = , 

27rT' ttT' ^ 



;i.3) 



where s is the entropy density, and T is the temperature. 

The finite coupling corrections were computed only for the shear viscosity. It was 
found in [20-23] that 



s Air 



^(n-™C(3)A-/= + ...). (1.4) 



In this paper, extending analysis of [20,21], we find 



2-ln2 375 ,, ^ , . ^ 

TnT = -—— + — C(3) A-3/2 + ■ ■ ■ , (1.5) 
/vr oZtc 



rj_ 

ttT 



(l-^C(3)A-^/^ + ...) . (1.6) 

The computations are quite technical, so we present only relevant steps and for 
the details refer the reader to previous work on the subject: [3,4,20,21,23]. In the 
next section we describe 0{a'^) near-extremal D3 brane geometry [24,25], primarily 
to set-up our notation. In section 3, following [20,23], we compute {rn, k} from the 
retarded correlation function of the stress energy tensor. In section 4, following [21,23], 
we compute the dispersion relation of a sound quasinormal mode up to third order in 
momentum, and confirm the value of rn obtained in section 3. 



^This is in fact universal to all four dimensional CFT gauge theories with a string theory dual. 
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2 Background geometry 



The background black brane geometry dual to a strongly coupled A/" = 4 SYM plasma 
at finite 't Hooft coupling was found in [24,25]. The ten dimensional background takes 
the form 



F5 = + *J^5, = —Advol 



(2.1) 



55 



where is a volume element of a round five-sphere, and 



dsl = g^audx^dx'' = ^e"(") i-fe^^'^Ue + dx^) + — ^e^^'^ . (2.2) 
u Au^ J 

Here f{u) = 1 — -u^, tq is the parameter of non-extremality of the black brane geometry, 
and we set the "AdS radius" L to one. To leading order irj^ 

7 = ^C(3)a'^ , (2.3) 

functions a, 6, c, v were found in [24, 25] 

a{u) =-157 (5^2 + hu^ - 3^6) , 
h{u) =15 7 {hu^ + 5m^ - 19m^) , 
c(ix) =0 , 

The Hawking temperature corresponding to the metric (12.1 p is 



(2.4) 



T = To (1 + 157) = - (1 + 157) . (2.5) 

TT 

These corrections were found assuming the only relevant term at order 7 is [28]. 
In [29] the full set of 7 corrections were computed including five-form terms, and there 
it was found that the black D3-brane solution receives corrections only from the 
term. In [17] it was further shown that the full spectrum of quasi-normal modes in this 
background is also unaffected by the five-form terms, and therefore we are justified in 
this paper in working just with the correction. 



*The string tension a', or more precisely a' /LP', is identified with A of tlie M — A SYM. 
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3 Relaxation time from the Kubo formula 



To obtain retarded correlation function of the boundary stress energy tensor, we study 
scalar perturbations of the background geometry (12. 2p (see [3,4]): 



It will be convenient to introduce a field ip{u, x), 

(p{u,x) = ^ h^y{u,x) , 

and use the Fourier decomposition 



(3.1) 



(3.2) 



ip{u,x) 



^—iujt+ik-x , ^ 



Finally, we introduce 



to 



k 



k 



2tiTq 2nTo 
3.1 The effective action 

The effective action to order 0(7) for ipk{u) takes form [20]: 



(3.3) 



(3.4) 



f d'k r\ 

=8^ y (2^ ^ 



+D ifkip^k + E ifilip^ + F ^lip'_ 



(3.5) 



The coefficients A, B, C, D, E, F are even functions of the momentum. They are given 
explicitly in appendix Rl 
Variation of Seff leads to 



iV2 

6S, - ' 



d^k 



[ du (EOM) + {Bi5ip^k + B2S^'_k) 
Jo 



87r2 J (2^)4 

where the coefficients of the boundary term are given by 



(3.6) 



i3i = -{Aif,)' + 2B^[ + C^u - 2{E^iy + F^l - {F^'J , (3.7) 
B2 = A^u + F^'^ + 2E^l , (3.8) 
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and EOM denotes the left hand side of the Euler-Lagrange equation 

Aipl + Cipi + 2D^k-^i2Bip', + dp, + FifD + ^iA^, + 2E^'l + F^',) = . (3.9) 

In order to have a well-defined variational principle, one has to add a generalized 
Gibbons-Hawking boundary term to the action (13.51) . As explained in [20], this should 
be done perturbatively in 7. Specifically, if we rewrite (13.91) in the form 

^l+Pi^'u+Po^k = 0{^), (3.10) 

with all 7-dependent terms exiled to the right, the generalized Gibbons-Hawking term 
fCgen rendering the variation of (13. 6p . takes form 

Kgen = -Aipk(p'_k - —v'kV'-k + Epi(p'i,ip'_k + 2EpoipkV'-k ■ (3-11) 

Notice that JCgen differs from the standard (supergravity) Gibbons- Hawking term JCstd- 

ICstd = -Aipk(p'_f, + /CiV9fcV9'_fc + fC2^k^-k , (3-12) 

where the coefficients /Ci oc 0(7) and IC2 oc 0(7°) are given in appendix Rl 
The bulk action (13. 5p can be rewritten in the form 



where 



A' C 

B = - Y^k^-k + B^W-k + -^VkV-k - E'^l^-k + ELpl^'_k - Eip'l^ip^i 

+ i^'pW-k - Y'^W-k ■ 



(3.14) 



Thus on-shell, it reduces to the sum of two boundary term: the horizon contribution 
( as M — s> 1 ) and the boundary contribution ( as m — ). In computing the two- 
point retarded correlation function of the boundary stress-energy tensor, the horizon 
contribution must be discarded [26]; the boundary contribution is divergent as u = 
e — s> and must be supplemented by the counterterm action [27]: 

^ct = - ^ ^ ^ d'x^ {^ + \P-^ [P'^Pki - P') In , (3.15) 
where jij is the metric induced at the u = e boundary, and 



P = f'P^J > P^^ = \ [R^, ' ^^7., ) • (3- 16) 



We can further rewrite (13.151) a: 



(3.17) 



where T is given in appendix |Al 

We would hke to emphasize that the counterterm action (I3.15P is computed assum- 
ing the standard Gibbons-Hawking term (13.121) . In particular, it takes into account the 
divergent as -u = e ^ boundary contribution coming from ( the supergravity part of 
) /C2 term in (I3.12p . Thus, even though at the level of the effective action (13. 5p the 
variational principle is well defined just with JCgen, in order to obtain finite retarded 
correlation functions with the counterterm action (13.150 . we need to supplement JCgen 
with the supergravity partO of )C2(pk'^-k- 



/C 



gen 



(3.18) 



Altogether, the total renormalized boundary action takes the form 

r d^k 



5, 



tot\ 



7r2 J {2TTy 



where 



•Fk 



(3.19) 



[B - A) (f[if_k + ]^{C - A + IC2 + T) (fk^_k - E'ifl^p.k + E^l(f'_^ 



F' 



WkV-k - yV'fcV'-fc + Epiip[ip'_i^ + 2Epoip[ip^k 



(3.20) 



'^Note that the conformal anomaly term in (|3.15|) does not contributed to the retarded correlation 
functions at order 0{w^,k^). 

^As observed in [20], the 0{j) parts of the coefScients in ([3ll|) . ([3?T2l) (and also (f3T7|) ) will not 
contribute in the limit u = e 0. 
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3.2 The solution for tpk 

We now turn to finding the solution to the equation of motion ( ]3.9p . or rather its 
equivalent to order 0(7), (13.101) : 



„ + 1 , k'^u^ + ra^ - fc^ 1 

"Pk + T2 = —.1 



3171m* + 3840k'u-^ + 2306n' 



- 600 j uip'k + ^ |^600tt)2 - 300fc2 + 50m + (3456/2^ - 2856w^)u^ + 768u^k^ (3.21) 

+ (-6560fc2 + 2136tt)2)M* + (-768fc* - 275)m^ + 3404fc2M'' + 225M^j ipk 
The incoming wave boundary condition for Lpk at the horizon implies [23] 

r, in:i(l — 157) 

ipkOc{l-u^) ^ — , u^l_. (3.22) 

Thus we represent the solution to (I3.2ip perturbatively in both 7 and (ro, fc^) as 



^k =(1 - ur""^ i + ^rv G^^ + + 0(tt)^ mfc^; 



(3.23) 



=g1^'°^ + 7Gf ^ + 0(f ) , ^ = 0,l,2, 
where 

Gl'^(ti) = G'i'^(ra,(?V) = A° Gl*^(Att),A2fc2;M), (3.24) 

and 



Explicitly, we find: 



lim (m) = 1 , lim GX\u) = lim ^(u) = . (3.25) 

It— >1- U^l— 11— >1_ 



Gr' = l. Gr = -|«=-|«^ + f , (3.26) 
Gr=0. Gr = f«« + ^«= + ^«'-f , (3.27) 
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G 



(2,1) 
k 



25 4 25 6 215\ / (u l\ 1, A /215, ^ 2445 



195 2 

— " + 



32 16 / 



369 25 



16 32 



In 2 



. 16 

^ u^^ ln(l + u) 



2445, „ 2885 215, 2. 605 /195, „ 195 

In 2 H In^ 2 u + In 2 

8 6 32 2 V 2 2 



2 40 3 
3 



+ 



135 25 



64 



In' 2 + 



1105 



43 



16 + 



, 43 25, 2. 369 

+ ln'2H 

2 64 16 



6 , f / 1375 ^ 25^6^ 25^4. j^^^^^^ 

1 e 16 ' 



In 2 U + ^ 



ra2 \ V 8 



16 



+ 



+ 



^21 1375 



195 



250 



6 



In 2 - 175m + —M^ 

2 3 



-u 



l^_^ln2]M4_g5^5 
2 16 



43 _ 25 
y ~ 16 



In 2 U' 



(3.28) 



3.3 Coupling constant correction to relaxation time 

Having found the solution for a gravitational perturbation, we can compute the corre- 
lation function Gxy,xy{^,<l) by applying the Minkowski AdS/CFT prescription [26] 

2K, 



Gl.>,g) = lim^ 



(3.29) 



Explicitly we find 



- — itv 



1 + I2O7 
905 



+ 



-q^ + - tt)^ln2 



+ 7 ( -120tt)2 In 2 + 25q2 + — ra^ 



(3.30) 



+ 0(ra^mq') +0(7^ 



where we used ( 12. 5 p to reintroduce the temperature, and denoted 



ro = 



q = 



(3.31) 

2nT ' 2-kT ^ ' 

In the hydrodynamic limit the retarded correlation function {u, q) takes form 



[4] 
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Comparing (I3.30p and (13.321) we conclude 

p = (i + 157 + o(f )) , - = i (i + 1207 + 0(f ; 

^ (3 33) 

^"^=^ + lr^ + °<^''- - = ^(l-l«7 + 0(7^)). 
Eq. (13.33P is our main result. 

4 Relaxation time from the sound pole 

The equation of motion for the sound quasinormal mode has been obtained in [21]. We 
represent it here in an equivalent form, to order 0(7): 

^"""'^ x(3m2 - (x2 + 2)q2) 
3tt)^ - (4x2(1 - x2)3/2 + 2tt)2(2x2 + l))q2 + x2(x2 + 2)q^ (4.1) 

^ (3m2 - (x2 + 2)q2)x2(l -x2)3/2 



~l~ 'Jsound[^sound\ ~l~ ^^(7 



2n 



where .^so^nd = Zgoundix), X = Vl — "^2 and (rD,q) are the dimensionless momenta 
reduced with respect to 27rTo, (13.41) . The order 0(7) the source term Jsound is given by 



'J sound ~ 



4x2(^3tt)2 - (x2 + 2)q2^ (l-x2)i/2 



^ sound Jsoundfl Zgound Jsound,! 

(4.2) 



Jsoundfl =256q2(l - x^)^(3tt)2 - (x^ + 2)q2) |^15ra^ - 30^2(5x2 + 2)q2 

+ (116x2 + 44 + 35x^)qM + ( 27tt)^(3171x^ + 4877 - 8648x2) 



+ 27tt)^(20291x2 - 12394 - 9790x^ + 3293x^)q2 - 9ro2(-24153x^ + 54172x2 
- 31188 + 51661x^ - 47492x^)q^ + (-322145x^ + 220694x* - 77416 

+ 89709x^° - 52806x^ + 147364x2)q6^ (1 - x^f^ , 

(4.3) 
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Jsound,! =648rD^(-5 - 59x2 ^ ^g^4^ _ 36ft,6(^_3939^4 _ ^gg ^ 7003^^ - 1259x2)q2 

+ 12tt)^q^(-40299x^ + 6242x^ - 360 + 159x^ + 40333x^) - 12tt)2q^(37927a;^ 

- 80 + 8900x2 _ 30499a;6 ^ ijogg^^iQ - 29972x^) + 4:x'^{x'^ + 2)(5811x^ 

+ 18043x6 - 14991x'^ - 12192x2 + 4004)q^ + |^-675x2tt)'^(9x^ - IGx^ + 5) 

+ 9x2tt)^q2(8190 + 87357x^ + 12989x2 - 111386x^) + q^((-79980x^ 

- 820917x^° + 1316595x^ - 284748x^ - 83700x2)tt)2 + 20736x2(1 - x2)ro^) 
+ q6(x2(131316x2 + 134283x^° + 24040 + 410574x^ - 430386x^ - 292777x^) 

+ (20736x6 + 20736x^ - 41472x2) ra^^) + (27648ro2x2 - 20736x6it)2 

- 6912x«ro2)q8 _ 768qiV(l - x2)(x2 + (1 - x2) V2 . 

(4.4) 

The incoming wave boundary condition at the horizon imphes the following perturba- 
tive expansion in the hydrodynamic limit [23] 



Z ^ -in)(l-157) ( ^(0) , (1) , g2 (2) ^/ 3 ^ 

^sound — I Zgound ^ '■^'^sound ^ ^ '^sound ^ ' ™1 

Si) _ (i) ^Oi'y'^') 
^ sound ~^sound,0 ^ l^sound,! ^ '-^W 7 ) 



(4.5) 



where 



and 



^iL.(^) - 4L(t^, ^) = A° 4L(A"), A2q2; X) , (4.6) 



iiSV ^^"""'^^^^ " ^ ' }^+ ^^^o^ridi^^ ^ }% ^'Sundi^) = • (4-7) 

The dispersion relation for the sound quasinormal modes is obtained by imposing a 
Dirichlet condition on Zgound the boundary: 

lim Zsoundix) = <^ m = w(q) . (4.8) 
4.1 Sound quasinormal spectrum to 0(q2) 

The sound quasinormal mode and and its spectrum to the leading and the first sub- 
leading order in the hydrodynamic approximation was found in [21,23]: 

(0) _ 3tt)2 + (x2 - 2)q2 _ 2mqx2 

soundfi 2j^2 _ 2q2 ' sound,Q <^^2 _ 2q2 ' \ ■ ) 
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sound,! 



(2404 + 446^2 - 4164x^ + 2006x^) 



16(3tt)2 - 2q2)2 

- 3ra2q2 (1588 + 183x2 - 2072x^ + 1003x^) + 45ro^ (5 - Ax^ + x^)^ , 

^'^ = 8q(3tt!^-2q2)2 ^^^^^^^ + ' ^^^0^' + 1734x6) 

- 3rt)2q2 (-9744 + 5035^^ - 2604x^ + 867x6) 

30roqx2 



- 36ro^ (594 - 264x2 + 43x^) 



2q2 - 3tt)2 

(4.10) 



tf (q) = - + ^7) + 0(q^ f ) ■ (4.11) 

4.2 Sound quasinormal spectrum to C(q^) 

Here, we extend the analysis of the previous section to the next order in the hydrody- 
namic approximation. Since our ultimate goal is to determine the dispersion relation 

(14.81) to order O(q^), it is sufficient to use the zeroth order dispersion relation, i.e., 

1 (2) 
to set ro = -^q. This drastically simplifies the hydrodynamic equations for z^J^^^q 

(2) 

and -Zgo^^^i- We find it more convenient to solve the resulting equations using the 
u = \f\ — x? variable. 
Explicitly we find 



^oL,o(-) =y (I + - dilog + ) + + ^ 



+ ^(m + 2)(1-m) 



(4.12) 



1 



,2 \ „,2 



^TJundM =u^C, + ( 2 + y ln(l - u')\ C, + ^Uu) 



-M+yln(l-^.2))x,(^) 



(4.13) 



where 



, ^ r , 4 + t2ln(l-t2N 
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T{t) ={ 10r(t + l)(9027r - 43808r + 58722r 

- 23100) ^dilog Q + - ^ 1^2(1 + + + l)t(10000 + (376440 

- 115500 In 2)t^ + (293610 In 2 - 1104128)t^ + (-219040 In 2 + 847972)t^ 
+ (45135 In 2 - 179418)t^) ln(l +t)- 18440 + (-20000 In 2 - 12568)t 
+ (-20000 ln2 + 1912)t=^ + (543780 - 752880 ln2 + 115500 In^ 2)^^ 
+ (1143490+ 115500 In^ 2 - 752880 In 2)t^ + (-293610 In^ 2 - 1403332 
+ 2208256 In 2)t^ + (-2629768 + 2208256 In 2 - 293610 In^ 2)f 
+ (219040 In^ 2 - 1695944 In 2 + 946434)t^ + (1729061 + 219040 In^ 2 

- 1695944 In 2)t^ + (-45135 In^ 2 + 358836 In 2 - 188378)t^ + (-340799 



+ 358836 In 2 - 45135 In^ 2)t 
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(4.16) 

and the integration constants Cj are tuned to satisfy the horizon boundary condition: 

lim ^oL,iM = 0. (4.17) 

The latter is achieved provided 

Ci = ^X2(l), C2 = -^Ji(l). (4.18) 

If we denote 

the Dirichlet boundary condition (14. 8 p will lead to the following dispersion relation for 
the sound quasinormal mode 

1 .2/1 105 \ J 3-21n2 



V3 V3 3 7 ^ \ 6v^ 

X (4.20) 
+ ^ (-2758 + I2zfl + 1705 In 2 h j + 0(q^ 7') . 

We were unable to evaluate (14.181) analytically; numerically, we find 

zfl = 264.7598406 . (4.21) 



13 



4.3 Relaxation time from the sound quasinormal spectrum 

Second order relativistic hydrodynamics of conformal fluids implies the following dis- 
persion relation for the sound mode [4] 

u = c,q- iTq^ + - f c^rn - ^ V' + ^(^^) ' (4-22) 



(4.23) 



cA 2 

where 

2r] 



3sT 

Comparing and we find 

c, = + ■ 7 + 0(7') , FT = i- (^1 + I2O7) + 0(7') , (4.24) 

in agreement with the conformal equation of state at order 0(7), as well as in agreement 
with the ratio ^ as given by (13.331) . Additionally, we compute 

TuT = + (2425 In 2 - 3358 + 124',] ) 7 + ^(t') • (4-25) 

ZTT lOTT V ' / 



A required agreement between (13.331) and (14.251) provides a prediction for 



(2) 



^1,0 



2425 



prediction ^ "^2 



which is in excellent agreement with the actually numerical result (I4.2ip . Thus we 
have a highly nontrivial check on our analysis. 
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A Coefficients of the effective action, the Gibbons-Hawking 
term, and the count erterm 



2u 



8f{u) + -fu^ I -600 + 25u^ + (-44tt)^ + Uk'^)u^ + 1760w^ - Uk\^ 



(A.l) 



8u 



24/(ti) + 77/^ I -1800 + 179u^ - (768fc^ + 768n)^)u^ + 5424u^ + 768fc^u 



.2^,5 



3387m^) 



(A.2) 



4^2/ 



8/(m)(3 + u^) + I 600 - 825u^ + (-104fe^ + 104tt)^)M^ - 16945ii 



+ (384fc^ + 872tv^)u^ + 34005ti'' - 280fe^M' - 16835m^ 



(A.3) 



D- "0 



16/(m)^ + 8uf{u)n)'^ + ^u^ ( eOOtt)^ + 250m + (-25fe^ + 25ro^)M^ 



8#/2 ^ 

+ {9UkW + 296rD^ + 296fc^ - 2120)u^ + (-2400tx)^ + 825k^)u^ + (1570 

- 592fc^ - 9UkW)u^ + (1007tt)^ - 1575fe^)^x^ + (296fc^ + 2220)^^ + 775^^^^ 



1920m' 
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/Ci = - 7 I -25 + (44ro^ - 44fc^)M + 160^^ + 44fc^M^ - 135m^ I . (A.7) 



F ^^2r^u'^f{u) {11 - 37u'^) . 

2 AAU'i\.. I 1«n.,2 I /l/ll.2„,3 



(A.4) 
(A.5) 
(A.6) 



4r^ 

/C2 = ^ 

M^ 



- 2 + 157^2 2u*^ - 8m^ + 5 



(A.8) 



15 



4 

r = - 



6 + (-2^2 + 2k'^)u - - 2fc\^ 

(A.9) 



2mV1 - 

+ 157^2 (3^^ - - 5 I I 3 + (ro^ + fc^)^ - 3^^ - fc^w^ 
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